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Abstract 


The  fluid  motion  about  an  elliptic  cylinder  im¬ 
pulsively  set  into  translational  and  rotational  mo¬ 
tion  is  obtained  by  the  method  of  matched  asymptotic 
expansions  for  small  time  and  large  Reynolds  number. 
The  constraint  of  the  perturbation  model  is  that  the 
boundary  layer  thickness  and  the  distance  of  travel 
are  of  the  same  asymptotic  order.  It  is  found  that 
pitch-up  motion  or  rotation  accompanying  translation 
at  an  angle  of  attack  is  indeed  capable  of  prevent¬ 
ing  the  early  formation  of  a  leading  edge  separation 
bubble.  Even  before  evident  in  the  streamline  pat¬ 
tern,  the  incipient  separation  bubble  is  accompanied 
by  a  characteristic  vorticity  signature  in  the  vi¬ 
cinity  of  the  leading  edge  Chat  is  quite  different 
from  that  with  rotation.  Further,  the  onset  of  an 
adverse  pressure  gradient  is  displaced  rearward  from 
its  location  for  pure  translation.  The  pre-Kutta 
condition  lift  evidently  arises  with  the  local  ac¬ 
celeration  that  is  a  consequence  of  the  displacement 
effect  of  the  growing  boundary  layer. 

I.  Introduction  'x 


The  dynamic  suppression  or  retardation  of  sepa¬ 
ration  from  an  airfoil  undergoing  rapid  pitching  mo¬ 
tion  enables  the  airfoil  to  achelve  an  angle  of  at¬ 
tack  that  may  be  much  larger  than  the  steady  flow 
stall  incidence  angle  and  yet  retain  an  intact 
boundary  layer.  This  is  concomitant  to  the  phenom¬ 
ena  of  dynamic  stall.  During  this  motion,  and  for 
unsteady  motion  in  general,  the  presence  of  reversed 
flow  adjacent  to  the  body  need  not  have  any  particu¬ 
lar  significance;  it  does  not  signal  insediate  sepa¬ 
ration  in  the  Prandtl  sense  of  gross  departure  of 
flow  from  the  surface  with  the  consequent  complete 
alteration  of  the  flow  field.*  Yet  its  presence 
seems  to  be  implicated  in  the  eventual  breakdown  of 
boundary  layer  flow,  whether  it  occurs  via  the 
bursting  of  a  leading  edge  separation  bubble,  per¬ 
haps  fed  by  a  thin  layer  of  reverse  flow  that  has 
advanced  upstream  along  the  surface  to  the  neighbor¬ 
hood  of  the  separation  bubble,  or  by  the  thickening 
of  this  layer  itself.2  It  was  for  insight  into  the 
delay  of  the  onset  of  flow  reversal  and  Its  spread 
that  the  present  work  was  undertaken. 

The  approach  taken  is  on*  that  is  recurrent  In 
the  literature;  namely,  the  development  of  an  asym¬ 
ptotic  solution  for  boundary  layer  growth  on  a  cir¬ 
cular  or  elliptic  cylinder  In  a  short  time  interval 
following  the  impulsive  onset  of  motion.  The  proce¬ 
dure  is  to  iterate  upon  a  basic  solution,  in  this 
case,  the  Rayleigh  solution  for  the  (locally)  flat, 
impulsively  started  plate.  By  continuing  in  this 
manner,  Coldstcin  and  Rosenhead^  were  able  to  extend 
Blaslus's  1908  work  to  a  third  approximation,  ob¬ 
taining  a  coordinate-type  expansion  In  small  tine 
for  the  boundary  layer  growth  on  a  circular  cylinder 
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(see  also  Batchelor*,  Stuart*).  Mutual  interaction, 
however,  between  the  growing  boundary  layer  and  the 
outer  flow  was  excluded;  the  pressure  distribution 
impressed  upon  the  boundary  layer,  for  Instance, 
remained  that  of  the  unperturbed  outer  flow.  Kang*, 
therefore  applied  the  technique  of  matched  inner  and 
outer  expansions  to  his  problem  of  the  Impulsively 
started  elliptic  cylinder.  Like  Goldstein  and 
Rosenhead,  he  was  able  to  trace  the  temporal  appear¬ 
ance  and  spatial  progression  of  reverse  flow  by 
solving  for  the  time  at  which  wall  shear  vanishes 
(if  it  does)  at  any  surface  location.  His  results, 
for  the  cylinder  moving  at  an  angle  of  attack,  show 
reversed  flow  advancing  upstream  from  the  vicinity 
of  the  rear  stagnation  point  and  eventually  merging 
with  a  slower-to-form  reverse  flow  region  near  the 
leading  edge  (a  separation  bubble).  For  the  present 
work,  we  adopt  the  method  of  Kang,  with  the  addition 
of  an  Impulsive  onset  of  rotation. 

II.  Perturbation  Model 


Initially  the  fluid  motion  is  everywhere  lrro- 
tatlonal  except  on  the  surface  of  the  cylinder  where 
viscous  interaction  requires  that  the  fluid  in  con¬ 
tact  with  the  surface  by  at  rest  relative  to  It  and 
where,  therefore,  a  surface  of  slip  or  vortex  sheet 
exists.  For  a  short  time  thereafter,  the  diffusion 
of  vorticity  is  the  dominant  effect  in  the  fluid 
close  to  the  surface.  It  is  this  short  time  inter¬ 
val  in  the  growth  of  the  boundary  layer  that  is  the 
concern  of  the  present  study. 

Posed  in  this  way,  it  is  perhaps  natural  to  con¬ 
sider  the  problem  as  a  singular  ceordinate-type  per¬ 
turbation  problem  for  small  values  of  the  independ¬ 
ent  variable  time  and  to  apply  matching  in  the  co¬ 
ordinate  normal  to  the  surface  to  overcome  the  dis¬ 
parity  in  the  inner  and  outer  characterizations. 
Although  there  is  nothing  wrong  with  such  an  ap¬ 
proach,  we  have  instead  chosen  to  formulate  the 
problem  as  a  parameter  perturbation  in  a  parameter 
c  which  tends  to  zero  as  time  tents  to  zero.  It 
represents,  in  the  domain  of  small  time,  a  rescaling 
of  the  geometric  time  scale  L/Lc  (characteristic 
length  and  speed)  analogous  to  the  R-*I  (for  Reynolds 
number  R)  rescaling  of  the  coordinate  normal  to  the 
surface  which  yields  the  familiar  boundary  layer 
equations. 


For  large  R,  another  source  o:  nonunlfomlty , 
and  a  practical  ccncem  in  evaluating  a  series  con¬ 
sisting  of  only  a  few  terms,  is  the  presence  of  ra¬ 
tios  of  the  temporal  perturbation  quantity  to  R. 

Kith  the  introduction  of  c ,  this  rmuniformitv  can 
be  expressed  as  c‘  i  o(c/R)  unlfomly  in  R  as  R  "•  **. 
It  becomes  convenient,  therfore,  to  establish  a  con¬ 
nection  between  the  limit  processes  of  R  *  *  and  c 
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The  alaplMt  connection  btdfftn  the  Halt  proc- 
essea  on  R  and  c  la  axpraaaed  by  R"1  «  ck  for  k  ».  0. 
For  a  apeclflc  value  of  k,  this  la  equivalent  to  an 
asynptotlc  order  relation  betveen  the  distance  d 
that  the  body  travels  in  time  t  at  a  speed  Uc  and 
the  quantity  4  ■  (vt)1*  <v  la  the  kineaatlc  vlscos- 
Ity)  which  characterlaea  the  boundary  layer  thick- 
nesa  for  unateady  fluid  motlca  about  a  body  atartlng 
froa  reat  at  t  •  0.  Thle  follows  from  letting  d/L  • 
0et/L  -  0(c)  and  writing  i2/l“  -  R“ld/L.  For  large 
k,  the  boundary  layer  assumptions  hold  to  a  good  ap¬ 
proximation,  that  la,  to  a  high  order  In  c.  For 
aaaller  k,  we  may  expect  an  earlier  departure  froa 
strict  boundary  layer  flow.  Let  us  take  the  model 
relation  R_1  «  c.  The  proportionality  constant  la 
then  the  local  Reynolds  nusber  based  on  the  length  d 
and  the  boundary  layer  thickness  is  of  the  sane 
asymptotic  order  as  the  distance  of  travel. 

111.  Formulation 

Let  us  normalize  the  equations  of  motion  with 
respect  to  the  length  L  taken  as  half  the  distance 
between  the  focli  of  the  elliptic  body,  the  density 
of  the  incompressible  fluid,  and  the  characteris¬ 
tic  velocity  combination  0C  ■  U#  +  (IL  of  the  speed 
of  the  distant  flow  Us,  incident  at  an  angle  of  at¬ 
tack  aQ,  and  the  rotational  velocity  flL  (fig.  1). 

For  convenience  we  shall  define  L/Ue  •  y  so  that 
Us/Uc  •  1  -  Y.  If  In  addition,  we  let  T  be  the 
characteristic  time,  representing  the  small  time 
scale  of  observation,  the  dinensionless  parameters 
UCT/L  and  R  are  found  to  completely  characterize  the 
fluid  motion  for  a  given  geometry  of  initial  and 
boundary  conditions.  In  accordance  with  the  previ¬ 
ous  discussion,  we  take 

t  -  UCT/L  -  B-2*-1  (1) 

where  S  la  a  constant  of  order  unity.  Then  the 
equations  of  motion  for  Incompressible,  two  dimen¬ 
sional  flow  in  a  body  fixed  frame  and  elliptic  cy¬ 
lindrical  coordinates  (C.n.z)  became 

it  h  J(  h'dt  h  d*J 


-el& 


xrev  - 

"1 


iv  r  f  id  ♦  3!  ♦  uw  ^  ^  1 

dt  *  I  h  df  h  d*l  h*  di|  h*  J 


We  shall  aasunc  that  the  outer  dependent  varia¬ 
bles  (lower  case)  have  expansions  of  the  form 

u(C,n,t;c)  “v  u0(c,n)  +  cut(C,n,t) 

+  c2u2(f,n,t)  +  .  .  . 
v  %  v0  +  cvj  +  e*v2  ♦  .  .  .  (5) 


Jkv  ^  Jhu  —  q 


where  u  *  uen  +  ve.  and  .the  dimensionless  scale  fac¬ 
tor  h  •  (sinhJ£  +  sin2n)’1.  The  last  terms  of  the 
momentum  equations  are  respectively  the  corlolls  and 
centrifugal  contributions  of  the  rotating  reference 
frame.  In  the  latter,  the  radial  position  vector 
has  been  expressed  as  r  ■  e-> h/3C  -  en3h/3n-  On  the 
surface  of  the  cylinder,  takes  to  be  the  coordinate 
surface  £  •  £fl,  u  *  v  ■  0.  As  •;  ■»  «  and  for  t  ■  0, 
the  solution  must  match  with  that  representing  a 
time  dependent  lnvlscld  flow. 


P  **•  P0  +  *P!  +  e*P2  +  •  •  • 

♦  %  *0  +  e*l  +  ‘S  ♦  •  •  • 

aa  C  ♦  0,  where  the  atreamfunctlon  satisfying  £q. 

(4)  Is  given  by  u  ■  h-I8<i/3£  and  v  •  -h-13p/3n. 

The  outer  flow,  being  initially  irrotatlonal  ex¬ 
cept  for  a  constant  background  vorticlty  of  magni¬ 
tude  2y  associated  with  the  rotating  frame,  remains 
strictly  Irrotatlonal  in  subsequent  approximations. 
Hence,  the  outer  fluid  motion  is  described  spatially 

by 

and  ’  2  <« 

»  0  for  n  »  0  (7) 

together  with  appropriate  boundary  conditions.  At 
t  ■  0,  “  0  for  n  >  0.  The  Bernoulli  forms  of 

Eqs.  (2)  and  (3), 

serve  to  determine  the  pressure  distribution. 

Th*  streanfunctlon  satisfying  (6)  and  describing 
the  initial  inviscid  flow  abouc  the  elliptic  cylin¬ 
der  is 

»j»  =  Sin  F£t) 

X  [tashll  -toshH.-tosHjJe^1  -l)J  . 

It  Is  normalized  to  zero  on  £  •-  and,  at  a  distance, 
in  the  original  (dimensional)  variables,  approaches 

-  (r-  5? )  v»(n*«.-At)  ♦  S(r*-r l) 

(let  r  -  Le^/2  and  (  ■*  •  while  L  •*  0  such  that  r  la 
fixed)  which  is  recognized  as  the  solution  for  the 
corresponding  flow  about  a  circular  cylinder  of  ra¬ 
dius  r0.  Expanding  $  for  small  c  ve  get 

f  =  O0) 

=  -0-r)e'*»inhU-l») 

*■  £  £toif>H  -  co»bl|,  -tell ij  (e ^  **-  l)J 


For  the  inner  expansions,  the  coordinate  normal 
to  the  surface  Is  rescaled  in  accordance  with  the 
expectation  that  a  length  scale  different  than  the 
geometric  length  L  is  more  appropriate  for  charac¬ 
terizing  local  variations.  This  is  accomplished  by 


2 


y 

applying  •  stretching  transformation  to  tha  coor- 
dlnata  {;  tha  raqulrad  atratch  balng  that  which  la 
aufflclaat  to  rastora  a  vlacoua  presence  In  tha 
Halt  t  +  0  (cf.  aqu.  (2)).  Define  tha  .atratchad 
variable  X  by 

C  -  C0  -  c  ?  (11) 

and  aaauaa  that  tha  inner  dependant  variablea 
(upper  case) ,  as  functions  of  X>  have  expansions 
of  the  font 

u(C.h.tit)  s.  uoU,n,t)  +  eUj^.n.t)  +  ... 

v  *  e{Vo  ♦  eVj  ♦  ...) 

p  n.  Pft  +  e?A  +  . . .  (12) 

♦  -v  c{f  +  t»  +  ...  ) 
o 

as  c  •*  0.  Further,  express  the  scale  factor  aa  a 
function  of  X  by  expanding  h  in  s  Taylor  series 
for  saall  c. 

a  h.(»t)(i*H h,(i) «••••] 


IV.  The  Solutisr.  for  R~X-  0(c) 

The  method  of  solution  is  formally  the  same  as 
that  of  Wang.  He  shall,  therefore,  merely  pres¬ 
ent  a  sketch  of  the  procedural  details.  Equations 
for  the  inner  tangential  T«loclty  and  outer  stream- 
function  perturbations  are  alternately  solved. 

Where  lacking,  boundary  conditions  and,  in  the  case 
of  the  inner  solution,  tangential  pressure  varia¬ 
tion  are  supplied  through  matching  conditions.  For 
this,  van  Dyke's  asymptotic  matching  principle  la 
adequate.  Step-by-step  we  may  also  check  and  con¬ 
firm  our  choice  of  an  asymptotic  sequence. 

The  zero-th  order  outer  solution  la  (cf.  equ. 

(10)) 

♦o  "  *o  <w> 

and  tha  corresponding  velocity  components  are 

11,  •  7-{-(i-r)e1,tosR(g*|#)sin(n-«.) 

1  ,  (21) 
4  JCf  ji*hH*teSH|e*‘* 
and  ’ 


Upon  substituting  the  expansions  for  u,  v,  and 
p  (Eq.  12)  into  the  governing  equations  (2)  and  (3) 
and  equating  like  powers  of  c,  we  get,  for  the 
first  inner  approximation. 


dU.  £U.  .  0 

dt  K  JV 

if.  =  0 
dX 


(13) 

(14) 


0»,a 


«*•  (»**♦) 

+  |s:«27(el(5  ?*-i)  J 


(22) 


The  solution  of  Eq.  (13)  for  the  zeroth  order 
of  the  inner  expansion  is 


U.=  U0(*.in)  er*  ;  '*xr,h 


(23) 


and,  for  the  second  approximation. 


au, £y, 

dt  hi  dV 


(15) 


.i  if, 

ho<)f 


if. 

dt 


p'illi  &) 

p  h,dijih«d^/ 


+  2  YU,  -  y'KsH, 


(16) 


Boundary  and  initial  conditions  on  UQ  and  U^  are 

Uo  ■  0  at  T  -  0,  t  >  0  (17) 

U  matches  outer  solution  as  (-*•,  t  »  0 
o  ^ 

and  for  t  ■  0,  (  >  0 

and  U^  *  0  at  C  -  0  and  t  •  0  (18) 

Uj  matches  outer  solution  as  (  -*  • 

The  streamfunctlcn  (zero  on  surface)  is  determined 
by  integrating  the  tangential  velocity  component 
in  the  direction  normal  to  the  surface: 


4*  a  j'uhil 

**  l  ,\  U9) 

•w  t[h,  J  J 


where  the  factor  u0(C0,r,)  has  been  identified  by 
matching  inner  and  outer  tangential  velocity  compo¬ 
nents.  U„  is,  of  course,  Rayleigh's  solution.  The 
streamfunctlon  in  this  approximation  is  obtained 
from  Eq.  (19). 

¥#  (24) 


The  surface  boundary  condition  for  Pj  is  found 
from  matching  the  inner  and  outer  streamfunction  ex¬ 
pansions,  (12)  and  (5),  with  pg  and  V0  given  by  Eqs. 
(20)and  (24).  This  yields  the  matching  condition 


U.(!.,*l)  (25) 


Pj  is  then  obtained  as  the  solution  of  the  Laplace 
equation  with  the  additional  condition  that 

♦  j  “  Pj  +  t*5’ (bounded  quantity)  as  (  ■»  ». 

Both  the  condition  on  the  surface  and  that  at  infin¬ 
ity  satisfy  the  spatial-tcrporal  separability  re¬ 
quirement  and  vanish  at  t  *  0.  Since  pj  ■  0  for  C 
■  C0.  the  solution  for  the  first  order  perturbation 
to  the  outer  streamfunctlcn  can  be  written  as 

P,  -  P,  +  C,  (26) 

with  p)  satisfying  (15)  and  expressed  as  a  Fourier 
series , 

1  /ft  v  2  mbi  J  W7) 
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whose  coefficients  say  be  evaluated  recursively 
froa  the  complete  elliptic  Integrals  of  the  first 
and  second  kind  vlth  the  parsaeter  depending  on  the 
eccentricity  of  the  elliptic  cylinder.  Conver¬ 
gence  of  the  series  Is  quite  rapid  for  fatter 
shapes  but,  as  expected,  auch  slower  for  eccentric¬ 
ities  approaching  unity. 

Tha  fluid  surface  upon  which  the  lnvlscld 
streaafunction  vanishes  to  first  order  no  longer 
coincides  with  the  surface  of  the  cylinder.  It 
has  been  displaced  by  the  growing  boundary  layer 
as  if  the  body  were  deforcing  In  tine  and  pres¬ 
enting  a  larger  obstacle  to  the  lnvlscld  fluid 
notion.  We  will  return  to  this  later  in  connec¬ 
tion  with  the  pressure  distribution. 

The  first  order  perturbation,  0^,  to  the 
Inner  tangential  component  of  velocity  Is  the  solu¬ 
tion  of  equ.  (15),  an  inhosogeneous  diffusion  eq¬ 
uation,  with  the  accompanying  Initial  and  boundary 
conditions.  Recall  equ.  (1-)  that.  In  the  first 
approximation,  the  pressure  was  constant  across  the 
boundary  layer.  Hence,  an  expression  for  3P  /3n  is 
found  by  matching  with  the  outer  expansion  for  this 
derivative  of  pressure.  Froa  equ.  (8),  this  gives 


dP. 

an 
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Performing  the  integrations  Indicated  in  equ.  (19), 
we  obtain 

t&(M>  ♦****( VI)  (31 

c-  =  *(»?-«)***  ♦  i  *') 

-l/TrF[q('  +  l**)Crf*  l}] 

-  Uef  f  U(2^3)er{5+r«?<r?,]j 


With  this  and  ocher  substitutions,  the  equation  for 
Uj  becomes 

35  ’H  <»> 

where 

EW  =  +&? e^5  ♦$«*(«*-») 

Foil  *  (  t*\ 

- 


uo  *  uo('o*n)  •  uo  *  iWn)/dn 

The  solution  prescribed  by  the  inhomogeneity 
has  the  form 


u,  =•  t  -r  t^zu,*}) 


When  the  boundary  condi  tions  y(0,n)  *  z(0,n)  are 
applied,  there  still  remains  an  unknown  function 
of  n  in  each  of  the  sol.utinns  y  and  z  which  is 
determined  from  catchim;  with  the  outer  expansion 
for  the  tangential  velocity  component.  The  results 
are  , 

-  i  -it-  if*  4l) 

eU  £?']1 

yfh 
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The  surface  condition  for  the  second  order 
perturbation  to  the  outer  streamfucction  is  easily 
obtained,  again  by  matching,  but  the  solution  Is 
complicated  by  the  presence  of  the  function  F(n) , 
Itself  expressed  as  a  series.  The  trend  is  to  in¬ 
creasingly  more  cumbersome  perturbations,  equations, 
and  solutions.  We,  therefore,  have  terminated  the 
development  at  this  stage,  though  it  turns  out  to 
have  been  perhaps  an  unfortunate  decision.  Since 
there  are  terms  of  the  inner  solution  having  no 
e2-match  (their  match  being  precisely  the  surface 
condition  on  C<2).in  the  composite  expansion,  the 
inner  solution  will  have  a  nonexponentially  decaying 
Influence  at  large  distances  from  the  body.  Never¬ 
theless,  this  influence  is  of  0(e2) 

A  uniformly  valid  composite  expansion  is  ob¬ 
tained  by  adding  the  inner  and  outer  streamfunctlon 
expansions  and  deleting  the  matching  conditions, 
which  occur  twice  in  the  sum.  Thus, 

*c  ■  ♦„  +  c<*,  +  +  e(*0  ♦  t?x)  -  **  (32> 

in  which 

♦U-t.fiffft),.  -  «(X-X.)h.u,U.(»i) 

consists  of  a  set  of  the  matching  conditions  and 
*c  is  expressed  uniformly  in  the  outer  variable  (. 
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Instantaneous  vorticity  contours,  vhers  chs 
vor deity  is  given  by  V2#  ,  are  overlaid  on  ths 
Instantaneous  streamlines'"  In  Figs.  2-5  for  an 
ellipse  of  eccentricity  2/3,  a  Reynolds  nuaber  of 
approximately  1000,  a  15°  initial  angle  'of  attack, 
and  at  a  dimensionless  tine  V  t/L  *  0.6.  y  “  0.4 
for  the  rotating  ellipse — a  qSite  high  pitch  rate 
0  L/U  "  2/3) — which  gives  an  Instantaneous  angle 
of  attack  of  v  28.7°  at  tine  0.6.  The  circled 
numbers  In  these  drawings  show  the  approximate 
location  of  local  surface  vortldty  extrema. 

We  have  found  that,  since  the  time  domain  of 
validity  shrinks  as  the  local  curvature  Is  in¬ 
creased,  the  only  appreciable  effect  of  Increasing 
the  eccentricity  Is  to  reduce  the  size  of  the  flow 
features;  when  drawn  scaled  to  the  same  tip  radius, 
the  patterns  are  essentially  Identical.  Further, 
as  time  advances,  flow  structures  merely  continue 
to  grow.  For  example,  the  trailing  edge  vortex 
shows  no  Inclination  to  depart  and,  in  trials  with 
a  non-rotating  circular  cylinder,  no  secondary  re¬ 
circulating  flow  structures  formed.  Obviously 
convective  Influences  are  not  strongly  represented 
In  a  diffusion  dominated  solution.  For  this  par¬ 
ticular  ellipse  and  Reynolds  number,  improbable 
diatortlon  in  the  flow  field  appeared  for  times 
greater  than  1.0  or  1.2.  Direct  comparison  with 
experiment  might  only  be  feasible  for  considerably 
smaller  times. 

Streamlines  and  vortldty  contours  for  the 
rotating  cylinder  (Figs.  4,5)  are  depicted  relative 
to  the  body  fixed  frame.  This  casts  a  slight 
shadow  on  the  interpretation  of  streamline  patterns. 
However,  the  interpretation  of  flow  structures, 
such  as  a  separation  bubble,  remain  fairly  reliable, 
especially  when  observed  in  a  continuum  of  rotation 
rates,  although  the  particular  form  may  vary.  This 
indicates  a  need  to  incorporate  all  available  de¬ 
scriptions,  including  the  pressure  and  vortldty 
fields,  into  interpretations.  Transformation  from 
the  body  frame  to  one  that  is  Instantaneously  co¬ 
incident  but  non-rotating, 

*  -  *.  -  i  (coah2C  -  sin2n)  , 

nr  d  l 

accomplishes  little  in  the  way  of  aiding  the  inter¬ 
pretation  of  streamllnea  but  it  shows  that  tbe 
pattern  of  vortldty  contours  is  invariant ,  tbe 
values  only  differing  by  2y.  (In  the  transformed 
system,  the  surface  of  the  rotating  cylinder  is  no 
longer  a  streamline;  streamlines  depart  from  that 
portion  of  the  surface  which  is  advancing  into  the 
fluid  and  are  drawn  to  that  portion  which  is  re¬ 
treating.) 

We  make  the  following  observations : 

1.  Without  rotation  (Fig.  3),  the  stress- 
line  pattern  is  a  deformation  of  the  early  sy=etric 
pattern  of  oscillating  flow  behind  a  circular  cy¬ 
linder.  At  the  Instant  depicted,  the  saddle  point 
has  been  pulled  apart  such  that  the  trailing  edge 
vortex  is  no  longer  bounded  by  a  closed  stagnation 
streamline  (a  piecewise  continuous  streamline  in 
the  fluid  whose  terminal  points  are  body  stagnation 
points);  the  vortex  is  free  to  depart  into  the 
wake.  Such  is  concomitant  to  the  establishment  of 
the  Kutta  condition  at  the  trailing  edge.  In 
contrast  to  this,  for  the  rotating  body  (Fig.  5), 
the  trailing  edge  vortex  is  confined  by  a  closed 
stagnation  streamline  (as  if  the  saddle  point  x^re 
pulled  apart  across  the  other  two  vertices)  er.d 
more  closely  resembles  a  separation  bubble.  Thus, 
the  establishment  of  the  Kutta  condition  la 


~) 


temporarily,  at  xeast,  suspended. 


2.  In  the  neighborhood  of  the  leading  edge 
(Figs.  2,4),  we  note  the  presence  of  very  high 
values  of  surface  vortldty  for  the  case  of  fixed 
angle  of  attack  as  compared  to  those  on  the  rotating 
ellipse.  This  is  one  Indication  of  the  fairly  sharp 
turn  the  flow  must  make  about  the  leading  edge  from 
the  stagnation  point  and  aupporta  the  view  that  the 
apparent  angle  of  attack  is  reduced  (Indeed  elimi¬ 
nated  in  the  case  drawn)  by  rotation.  Additional 
support  comes  from  the  presence  of  the  cell  of 
positive  vortldty  lying  just  off  the  leading  edge 
which  is  completely  absent  in  the  case  of  rotation. 
Its  presence  may  be  taken  as  a  precursor  or  sign 
of  an  incipient  leading  edge  aeparation  bubble  and 
conversely  for  its  absence.  (A  aeparation  bubble 
is  not  evident  from  the  streamline  pattern  for  such 
a  low  angle  of  attack.  At  v  28  ,  the  instantaneous 
angle  of  attack  of  the  rotating  body,  a  bubble  Is 
well  developed.) 


V.  Pressure  Distribution  and  Lift 

The  zeroth  order  surface  pressure  distribution 
is  first  that  of  the  outer  lnviscid  flow  impressed 
across  the  bl  (P„(n)  "  P0(C  ,n),  ct.  Eq.  (14)).  In¬ 
tegrating  Eqs.  (8,9)  along  a  path  running  from  a 
point  at  infinity  to  an  arbitrary  point  (C0,n)  on 
the  surface  of  the  ellipse  and  expressing  the 
result  as  a  pressure  coefficient,  we  get (for  y  i  1, 
i.e.,  Ug  i  0) 

Cr. «  *3.  =  i  -  [  -  j(  -  s,vn) 

.(if.;.,:).  £1^1 


(33) 


where  ^  is  a  velocity  potential  defined  by  u. 
for  the  time  dependent  first  order  (outer)  velocity, 

♦,  -  t  r(i-y) 


yxr  it-i 


in  which  the  first  term  stems  from  the  changing 
angle  of  incidence  of  the  distant  stream  and  the 
second  represents  the  displacement  affect  of  the 
growing  bl.  The  pressure  singularity  at  t  »  0  is 
traceable  to  the  impulsive  onset  of  motion:  The 
initial  rate  of  growth  of  the  BL,  estimated  by  the 
initial  pseed  at  which  vorticity  penetrates  into 
the  fluid  with  penetration  depth  taken  in  the  usual 
sense  as  that  depth  at  which  the  vorticity  has  a 
value  1/e  tines  its  surface  value,  is  singular  at 
t  “  0,  so  also  is  the  outer  inviscid  flow  response. 

In  order  to  determine  the  first  order  pressure 
perturbation,  inner  or  outer,  it  would  be  necessary 
to  know  the  second  order  perturbation  to  the  outer 
strcamfunctlon,  (>2 »  or  at  least  dvj/Bf  calculated 
at  (0.  This  should  not  be  too  difficult  to  calculate 
but  has  not  been  included  here. 

Lift  is  determined  from  the  integral 

SIT 

(cr(<l)-cr(l«'l)]  (34) 
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in  which  a(t)  is  the  instantaneous  angle  of  attack 
and  the  projection  of  surface  area  onto  the  chord 
is  given  by  the  ratio  of  |dx(o)|  ■  cosh£0  sinndn 
to  hdn.  This  is  easily  integrated  to  give 

tL- 

with 

A,  *  i\ 
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the  only  remaining  Fourier  coefficient. 

The  pressure  (C )  distributions  for  the 
rotating  and  non-rotating  cylinders  are  drawn  for 
comparison  in  Fig.  6  and  the  associated  values  of 
lift  for  various  times  are  presented  in  Table  1. 

For  the  rotating  cylinder,  we  notice  the  extended 
region  of  favorable  pressure  gradient  on  the  upper 
surface  which  qualitatively  resembles  the  distri¬ 
bution  on  a  nonrotating  elliptic  cylinder  at  a  very 
small  or  zero  angle  of  attack.  The  pressure  drop 
near  the  trailing  edge,  also  on  the  upper  surface, 
occurs  under  the  influence  of  the  trailing  edge 
vortex.  In  contrast,  the  non-rotating  body  ex¬ 
hibits  the  usual  extended  region  of  adverse  gradient 
and  diagonal  symmetry.  Lift  is  seen  to  be  sub¬ 
stantially  higher  for  the  rotating  ellipse  than  for 
that  at  a  fixed  angle  of  attack.  The  t-1/2  decay 
is,  in  this  approximation,  the  vestige  of  the 
impulsive  start. 

The  dimensionless  circulation  Inferred  from 
the  presence  of  lift  is  (^/(l  -  y) .  At  the  time  of 
this  writing,  we  have  not  succeeded  in  spatially 
Integrating  the  vortlclty  to  confirm  the  values  of 
circulation.  In  part,  the  difficulty  lies  in  the 
non-exponential  decay  of  vorticity  mentioned  pre¬ 
viously. 

VI.  Conclusion 

Within  the  limitations  of  the  solution,  we 
have  seen  here  the  addition  of  pitching  motion  or 
rotation  suppresses  the  early  formation  of  leading 
edge  separation  bubbles  by  creating  a  situation  in 
which  che  effective  angle  of  attack  is  much  reduced. 
We  have  also  found  the  Initial  lift  to  be  enhanced 
and  che  establishment  of  the  Kutta  condition  at  the 
trailing  edge  to  have  been  put  in  abeyance  by 
rotation. 
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Fig.  1.  Physical  configuration. 


Fig.  2.  Leading  edge  of  elliptic  cylinder, 
translation  to  left  at  15°  angle 
of  attack  at  time  Uct/L  •  0.6. 
Streamlines  (dashed)  and  vortlclty 
contours  (solid,  labeled)  are  over¬ 
laid.  Circled  values  indicate 
approximate  location  of  vortlclty 
extrema. 
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3.  Trailing  edge  of  elliptic  cylinder 
translating  at  13°  angle  of  attack 
at  time  0.6  (cf.  Fig.  2). 


.180  .127  .104  .090  .030 

4.104  3.091  3.704  3.510  3.301 


Trailing  edge  of  elliptic  cylinder 
rotating  and  translating  at  an 
instantaneous  angle  of  attack  of 
28.7°  (pitch-up  froa  II’),  for 
tine  0.6  and  rotation  rate  2/3 
(cf.  Fig.  2). 


Table  1.  Calculated  lift  versus  tine 


Chord  Position 


Leading  edge  of  elliptic  cylinder 
rotating  clockvlse  and  translating 
at  an  instantaneous  cnglc  of  attack 
23.7°,  tine  0.6,  and  rotation  rate 
0L/U-  -  2/3.  (cf.  Fig.  2) 


Fig.  6.  C  distribution  or.  elliptic  cylincci 
Dashed  curve  is  for  cylinder  trans¬ 
lating  at  13  angle  of  attack:  soli< 
curve  is  for  cylinder  translating  ai 
rotating  (pitch-up)  fret:  an  initial 
angle  of  attack  of  !i  .  Time  0.1. 
Upper  and  lower  surface,  are  dcr.otc< 
by  U  and  L. 


